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Abstract. Super-Kamiokande collaboration assumes that the direction of every 
observed lepton coincides with the incoming direction of the incident neutrino, which 
is the fundamental basement throughout all their analysis on neutrino oscillation. We 
examine whether this assumption to explain the experimental results on neutrino 
oscillation is theoretically acceptable. Treating every physical process concerned 
stochastically, we have examined if this assumption just cited is acceptable. As the 
result of it, we have shown that this assumption does not hold even if statistically. 
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1. Introduction: The motivation of the paper 

According to the results obtained from the Super-Kamiokande Experiments on 
atmospheric neutrinos, oscillation phenomena have been found between two neutrinos, 
v^&nd v T . Published reports on the confirmation to the oscillation between the neutrinos, 
z/^and u T , and the history forgoing to these experiments will be critically reviewed and 
details are in the following: 

(1) During 1980's Kamiokande and 1MB observed smaller atmospheric neutrino flux 
ratio z/^/Vg than the predicated value pQ. 

(2) Kamiokande found anomaly in the zenith angle distribution j2] 

(3) Super-Kamiokande found v lL -v T oscillation [3] and 

Soudan2 and MACRO confirm the Super-Kamiokande result pE]. 

(4) K2K, the first accelerator-based long baseline experiment, confirmed atmospheric 
neutrino oscillation [5]. 

(5) MINOS's precision measurement gives the consistent results with Super- 
Kamiokande ones [6]. 

It is well known that Super-Kamiokande Collaboration examined all possible types 
of the neutrino events, such as, say, Sub-GeV e-like, Multi-GeV e-like, Sub-GeV //-like, 
Multi-GeV /i-like, Multi-ring Sub-GeV //-like, Multi-ring Multi-GeV //-like, PC, Upward 
Stopping Muon Events and Upward Through Going Muon Events, in other words, all 
possible interactions by neutrino, such as, elastic and quasielastic scattering, single- 
meson production and deep scattering are considered. Furthermore, all topologically 
different types of neutrino events lead the unified numerical oscillation parameters, say, 
Am 2 = 2.4 x 10" 3 eV 2 and sin 2 26 = 1.0 [12]. 

Taking into account all factors mentioned above, it is natural that the majority 
should believe the finding of the \l — r neutrino oscillation by Super-Kamiokande 
Collaboration. 

However, it should be emphasized strongly that the Super-Kamiokande 
Collaboration put the fundamental assumption in the analysis of the atmospheric 
neutrino oscillation which is never self-evident and should be carefully examined. This 
assumption is that the direction of the incident neutrino is the same as that of emitted 
lepton(s). 

In order to avoid any misunderstanding toward the SK assumption on the direction 
we reproduce this assumption from the original SK papers and their related papers in 
italic. 

Kajita and Totsuka state: 

"However, the direction of the neutrino must be estimated from the 
reconstructed direction of the products of the neutrino interaction. In water 
Cheren-kov detectors, the direction of an observed lepton is assumed to be the 
direction of the neutrino. Fig. 11 and Fig. 12 show the estimated correlation 
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angle between neutrinos and leptons as a function of lepton momentum. At 
energies below 400 MeV/c, the lepton direction has little correlation with the 
neutrino direction. The correlation angle becomes smaller with increasing 
lepton momentum. Therefore, the zenith angle dependence of the flux as a 
consequence of neutrino oscillation is largely washed out below 400 MeV/c 
lepton momentum. With increasing momentum, the effect can be seen more 
clearly.' 1 [7J 

On the other hand, Ishitsuka states in his Ph.D thesis which is exclusively devoted into 
the L/E analysis of the atmospheric neutrino from Super-Kamiokande experiment as 
follows: 

" 8.4 Reconstruction of L v 

Flight length of neutrino is determined from the neutrino incident zenith an- 
gle, although the energy and the flavor are also involved. First, the direction of 
neutrino is estimated for each sample by a different way. Then, the neutrino 
flight lenght is calclulated from the zenith angle of the reconstructed direction. 

8.4-1 Reconstruction of Neutrino Direction 
FC Single-ring Sample 

The direction of neutrino for FC single-ring sample is simply assumed to be 
the same as the reconstructed direction of muon. Zenith angle of neutrino is 
reconstructed as follows: 

cos & r u ec = cos 6 M (8.17) 

, where cos6™ c and cos9 M are cosine of the reconstructed zenith angle of 
neutrino and muon, respectively. 17 [8\ 

Furthermore, Jung, Kajita et al. state: 

"At neutrino energies of more than a few hundred MeV, the direction of 
the reconstructed lepton approximately represents the direction of the original 
neutrino. Hence, for detectors near direction of the lepton. Any effects, such 
as neutrino oscillations, that are a function of the neutrino flight distance will 
be manifest in the lepton zenith angle distributions." [9] 

Hereafter, we call the fundamental assumption by the Super-Kamiokande Experi- 
ment simply as the SK assumption on the direction. 

In section 2, we discuss treatment on the cross section for Quasi-Elastic 
Scattering(QEL) stochastically which play an particularly important role in the analysis 
of Fully Contained Events. In section 3, we treat the effect of azimuthal angle of 
QEL over the zenith angle of the neutrino events which could not be treated in the 
Detector Simulation carried by the Super-Kaomiokande Collaboration. In section 4, we 
treat the zenith angle distributions of the neutrino events for the different directions of 
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the incident neutrinos stochastically, taking into account of the corresponding incident 
neutrino energy spectrum. In section 5, we discuss the correlation diagram between the 
zenith angles of the incident neutrinos and those of the produced muons, and show that 
the SK assumption on the direction does not hold even if statistically. As the result of 
it, for example, upward neutrino zenith angle distribution in the absence of neutrino 
oscillation does not reproduce the corresponding one of the produced muons. In the 
conclusion, we state that we will give the results of re-analysis of the L/E distribution 
obtained by Super-Kamiokande experiment under the situation that the SK assumption 
on the direction does not hold in the subsequent paper. 



2. Cross Sections of Quasi Elastic Scattering in the Neutrino Reaction and 
the Scattering Angle of Charged Leptons. 

In order to examine the validity of the SK assumption on the direction, we consider the 
following quasielastic scattering(QEL): 
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because the QEL is the most dominant process in Fully Contained Events and Partially 
Contained Events which occur in the detector. Also, we adopt the QEL events among 
Fully Contained Events as the object for the analysis, because they are simple and solid 
events in which the energies of the produced muons (electrons) could be uniquely de- 
termined, compared with any other type of the neutrino event. The differential cross 
section for QEL is given as follows |10j . 
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The signs + and — refer to io e ) and v^i e \ for charged current (c.c.) interactions, 
respectively. The Q 2 denotes the four momentum transfer between the incident neutrino 
and the charged lepton. Details of other symbols are given in |10j . 

The relation among Q 2 , E V ( V \, the energy of the incident neutrino, Eg, the energy 
of the emitted charged lepton (muon or electron or their anti-particles) and 6 S , the 
scattering angle of the emitted lepton, is given as 

Q 2 = 2E u{D) E e (l-cos6 s ). (3) 

Also, the energy of the emitted lepton is given by 

Q 2 

Ei = E v(p) - — . (4) 

Now, let us examine the magnitude of the scattering angle of the emitted lepton 
in a quantitative way, as this plays a decisive role in determining the accuracy of the 
direction of the incident neutrino, which is directly related to the reliability of the zenith 
angle distribution of both Fully Contained Events and Partially Contained Events in the 
Super-Kamiokande Experiment. 

By using Eqs. (J2J) to (jl]), we obtain the distribution function for the scattering 
angle of the emitted leptons and the related quantities by a Monte Carlo method. The 
procedure for determining the scattering angle for a given energy of the incident neutrino 
is described in the Appendix A. Fig. [1] shows this relation for muon, from which we can 
easily understand that the scattering angle 6 S of the emitted lepton ( muon here ) cannot 
be neglected. For a quantitative examination of the scattering angle, we construct the 
distribution function for 9 S of the emitted lepton from Eqs. ([2]) to (j4]) by using a Monte 
Carlo method. 

Fig. [2] gives the distribution function for 6 S of the muon produced in the muon 
neutrino interaction. It can be seen that the muons produced from lower energy 
neutrinos are scattered over wider angles and that a considerable part of them are 
scattered even in backward directions. Similar results are obtained for anti-muon 
neutrinos, electron neutrinos and anti-electron neutrinos. 

Also, in a similar manner, we obtain not only the distribution function for the 
scattering angle of the charged leptons, but also their average values < 9 S > and their 
standard deviations a s . Table 1 shows them for muon neutrinos, anti-muon neutrinos, 
electron neutrinos and anti-electron neutrinos. From the Table 1, it seems that the 
scattering angles could not be neglected. However, Super-Kamiokande Collaboration 
approximate them as zero, as cited in three passages (in Italic) mentioned above 
[7J,[8j,[9j. Or they may claim that such the assumption surely holds statistically as 
whole (see the Conclusion of our paper, relating to the latter surmise). 

3. Influence of Azimuthal Angle of Quasi Elastic Scattering over the Zenith 
Angle of both Fully Contained Events and Partially Contained Events 



In the present section, we examine the effect of the azimuthal angles, 0, of the emitted 
leptons over their own zenith angles, ^(/j)), for given zenith angles of the incident 
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Figure 2. Distribution functions for the scattering angle of the muon for muon- 
neutrino with incident energies, 0.5 , 1.0 and 2 GeV. Each curve is obtained by the 
Monte Carlo method (one million sampling per each curve). 



neutrinos, 6 U ^, which could not be considered in the detector simulation carried by 
the Super-Kamiokande Collaboration, [j] 

For three typical cases (vertical, horizontal and diagonal), Fig. 3 gives a schematic 
representation of the relationship between, 9 v {p), the zenith angle of the incident 
neutrino, and ( 6 S , <fi), a pair of scattering angle of the emitted lepton and its azimuthal 
angle. 

From Fig. 3(a), it can been seen that the zenith angle of the emitted lepton 
is not influenced by its <fi in the vertical incidence of the neutrinos (9 u (u) = 0°), as it 
must be. From Fig. 3(b), however, it is obvious that the influence of <fi of the emitted 

| Throughout this paper, we measure the zenith angles of the emitted leptons from the upward vertical 
direction of the incident neutrino. Consequently, notice that the sign of our direction is oposite to that 
of the Super-Kamiokande Experiment ( our cos#„( P ) = - cos#„( P ) in SK) 
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Table 1. The average values < 6 S > for scattering angle of the emitted charged leptons 
and their standard deviations a s for various primary neutrino energies E v ^y 



E U{P ) (GeV) 


angle 
(degree) 




"MA) 
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<8 S > 
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89.74 


67A7 
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0.5 


<9 S > 
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32.82 


1 


<e s > 


48.44 


36.00 


48.42 


36.01 
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32.06 


27.05 


2 


<e s > 


25.84 
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25.84 


20.20 


OS 


21.40 


17.04 


21.40 


17.04 
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<e s > 


8.84 


7.87 


8.84 


7.87 


OS 


8.01 


7.33 


8.01 


7.33 


10 


<e s > 


4.14 


3.82 


4.14 


3.82 


OS 


3.71 


3.22 


3.71 


3.22 


100 


<e s > 


0.38 


0.39 


0.38 


0.39 


OS 


0.23 


0.24 


0.23 


0.24 



leptons on their own zenith angle is the strongest in the case of horizontal incidence of 
the neutrino {0 v (p) = 90°). Namely, one half of the emitted leptons are recognized as 
upward going, while the other half is classified as downward going ones. The diagonal 
case ( 6 u (u) — 43°) is intermediate between the vertical and the horizontal. In the 
following, we examine the cases for vertical, horizontal and diagonal incidence of the 
neutrino with different energies, say, E v (p\ = 0.5 GeV, E v r v \ = 1 GeV and E v {p\ = 5 
GeV, as the typical cases. 

3.1. Dependence of the spreads of the zenith angle for the emitted leptons on the 
energies of the emitted leptons for different incident directions of the neutrinos with 
different energies 

The detailed procedure for the Monte Carlo simulation is described in the Appendix A. 
We give the scatter plots between the fractional energies of the emitted muons and their 
zenith angle for a definite zenith angles of the incident neutrino with different energies 
in Figs. H] to [6j In Fig. HI we give the scatter plots for vertically incident neutrino with 
different energies 0.5 GeV, 1 GeV and 5 GeV . In this case, the relations between the 
emitted energies of the muons and their zenith angles are unique, which comes from the 
definition of the zenith angle of the emitted lepton. However, the densities (frequencies 
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Figure 3. Schematic view of the zenith angles of the charged muons for different 
zenith angles of the incident neutrinos, focusing on their azimuthal angles. 

of the event number) along each curves are different in position to position and depend 
on the energies of the incident neutrinos. Generally speaking, densities along the curves 
become greater toward cos O^p) — 1. In this case, cos is never influenced by the 
azimuthal angel in the scattering by the definition l§|. 

On the contrast, it is shown in Figure 5 that the horizontally incident neutrinos 
give the widest zenith angles distribution for the emitted energies of muons due to the 
effect of their azimuthal angles. The more lower incident neutrino energies, the more 
wider spread of the emitted leptons. As easily understood from Figure 6, the diagonally 
incident neutrinos give the intermediate zenith angle distributions for the emitted muons 
between those for vertically incident neutrinos and those for horizontally neutrinos. 

3.2. Zenith angle distribution of the emitted lepton for the different incidence of the 
neutrinos with different energies 

In Figures 7 to 9, we express Figs. 4 to 6 in a different way. We sum up muon events 
with different emitted energies for given zenith angles. As the result of it, we obtain 
frequency distribution of the neutrino events as a function of cosO^ for different incident 
directions and different incident energies of neutrinos. 

§ The zenith angles of the particles concerned are measured from the vertical direction. 
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Figure 4. The scatter plots between the fractional energies of the produced muons 
and their zenith angles for vertically incident muon neutrinos with 0.5 GeV, 1 GeV 
and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Figure 5. The scatter plots between the fractional energies of the produced muons 
and their zenith angles for horizontally incident muon neutrinos with 0.5 GeV, 1 GeV 
and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Figure 6. The scatter plots between the fractional energies of the produced muons 
and their zenith angles for diagonally incident muon neutrinos with 0.5 GeV, 1 GeV 
and 5 GeV, respectively. The sampling number is 1000 for each case. 
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Figure 7. Zenith angle distribution of the muon for the vertically incident muon 
neutrino with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 
10000 for each case. SK stands for the corresponding ones under the SK assumption. 
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Figure 8. Zenith angle distribution of the muon for the horizontally incident muon 
neutrino with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 10000 
for each case. SK stands for the corresponding ones under the SK assumption. 



mu 

CO! 


Dn r 

)6 V = 


eut 
0.7 


'ino 
31(6 


v=4 


3°) 










a\ 


E v = 
'era 


0.5 

ge= 


3eV 
0.1 


39 






S 


K 






s 


d.= 


0.5E 


6 










































































-w-r- 















-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 
COS8.. 



(b) 



mu 

CO. 


Dn r 
S8 V = 


eut 
0.7 


nno 
31(6 


v=4 


3°) 










a\ 


E v = 
'era 


1G 

ge= 


3V 

0.4 


32 






S 


K 






s 


d.= 


0.4* 


53 









































































-0.8-0.6-0.4-0.2 0.2 0.4 0.6 O.t 
COS8„ 



mupn neutrino 

0.7:i1(^ v =4p 
SGbV 



CO£S6„H' 



average 



0.715 



d.= 



0.103 



-0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 
COS6 u 



Figure 9. Zenith angle distribution of the muon for the diagonally incident muon 
neutrino with 0.5 GeV, 1 GeV and 5 GeV, respectively. The sampling number is 10000 
for each case. SK stands for the corresponding ones under the SK assumption. 
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neutrinos with the vertical, horizontal and diagonal directions, respectively (see 
Figure 3). The overall neutrino spectrum at Kamioka site is taken into account. The 
sampling number is 10000 for each case. SK stand for the corresponding ones under 
the SK assumption. 



In Figures 7(a) to 7(c), we give the zenith angle distributions of the emitted muons 
for the case of vertically incident neutrinos with different energies, say, E u = 0.5 GeV, 
E v = 1 GeV and E v = 5 GeV. 

Comparing the case for 0.5 GeV with that for 5 GeV, we understand the big contrast 
between them as for the zenith angle distribution. The scattering angle of the emitted 
muon for 5 GeV neutrino is relatively small (See, Table 1) so that the emitted muons 
keep roughly the same direction as their original neutrinos. In this case, the effect of 
their azimuthal angle on the zenith angle is also small. However, in the case for 0.5 
GeV which is the dominant energy for Fully Contained Events and Partially Contained 
Events in the Superkamiokande, there is even a possibility for the emitted muon to be 
emitted in the backward direction due to the large angle scattering, the effect of which 
is enhanced by their azimuthal angle. 

The most frequent occurrence in the backward direction of the emitted muon 
appears in the horizontally incident neutrino as shown in Figs. 8(a) to 8(c). In this 
case, the zenith angle distribution of the emitted muon should be symmetrical to the 
horizontal direction. Comparing the case for 5 GeV with those both for ~0.5 GeV and 
for ~1 GeV, even 1 GeV incident neutrinos lose almost the original sense of the incidence 
if we measure it by the zenith angle of the emitted muon. Figs. 9(a) to 9(c) for the 
diagonally incident neutrino tell us that the situation for diagonal cases lies between the 
case for the vertically incident neutrino and that for horizontally incident one. SK in 
the figures denotes the SK assumption on the direction of incident neutrinos. From the 
Figures 7(a) to 9(c), it seems to be clear that the scattering angles of emitted muons 
influence their zenith angles, which is enhanced by their azimuthal angles, particularly 
more inclined directions of the incident neutrinos. 
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4. Zenith Angle Distribution of Fully Contained Events and Partially 
Contained Events for a Given Zenith Angle of the Incident Neutrino, 
Taking Their Energy Spectrum into Account 

In the previous sections, we discuss the relation between the zenith angle distribution of 
the incident neutrino with a single energy and that of the emited muons produced by the 
neutrino for the different incident direction. In order to apply our motivation around 
the uncertainty of the SK assumption on the direction for Fully Contained Events and 
Partially Contained Events, we must consider the effect of the energy spectrum of the 
incident neutrino. The Monte Carlo simulation procedures for this purpose are given in 
the Appendix B. 

In Fig. [TOj we give the zenith angle distributions of the sum of and \l~ 

for a given zenith angle of Dp, and z/^, taking into account different primary neutrino 
energy spectra for differnt directions at Kamioka site. SK in the figures denotes the 
SK assumption on the direction. From the figures, it seems to be clear that the 
SK assumption on the direction does not hold. Namely, we could conclude that the 
scattering angle of the emitted muons acompanied by their azimuthal angles influence 
their zenith angle distribution for the all directions. 

5. Correlation between the Zenth Angle Distribution of the Incident 
Neutrinos and that of the Emitted Leptons 

Now, we extend the results for the definite zenith angle obtained in the previous sections 
to the case in which we consider the zenith angle distribution of the incident neutrinos 
totally. 

Here, we examine the real correlation between cos^ and cos 9^, by performing the 
exact Monte Carlo simulation. 

The detail for the simulation procedure is given in the Appendix C. 

In order to obtain the zenith angle distribution of the emitted leptons for that of 
the incident neutrinos, we divide the cosine of the zenith angle of the incident neutrino 
into twenty regular intervals from cos9 u = to cos9 u = 1. For the given interval of 
cos8 u , we carry out the exact Monte Carlo simulation, and obtain the cosine of the 
zenith angle of the emitted leptons. 

Thus, for each interval of cos9 u , we obtain the corresponding zenith angle 
distribution of the emitted leptons. Then, we sum up these corresponding ones over 
all zenith angles of the incident neutrinos and we finally obtain the relation between the 
zenith angle distribution for the incident neutrinos and that for the emitted leptons. 

In a similar manner, we could obtain between cos 6^ and cos 9p for anti-neutrinos. 
The situation for anti- 
neutrinos is essentially the same as that for neutrinos. 

In Fig. [11] we classsify the correlation between cos 9 V and cos 9^ according to the 
different energy range of the incident muon neutrinos. The cos 9^ distribution along 
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Figure 11. Correlation diagram between cos 9 V and cos 9^ for different neutrino energy 
regions. 



cos 6^ = 1, cos9 u = and cos 6^ = 0.73 in Figure 11 correspond to Figure 10(a) 
(vertical), Figure 10(b) (horizontal) and Figure 10(c) (diagonal), respectively. Looking 
the cos^ distribution for the fixed cos 9 V in the Figure 11, it is well understood that 
the cos9 fjL distribution around cos6^ ~ cos 6^ spreads wider as cos 6^ decreases. This is 
due to the effect of scattering angle enhanced by the azimuthal angle (see Figures 3 and 
15, also). 

In Fig. [12} we classify the correlation between cos9 v and cos # M according to the 
different energy range of E^. It is clear from Figure 12 that (a)the SK assumption on 
the direction (cos 9^ « cos 6^) does hold surely > 5 GeV, (b)this assumption does 
hold rougly > 2 GeV. However, as E^ decreases, this relation becomes incorrect. In 
the energy range of 0.5 < E^ < 1 GeV where neutrino events in the Super-Kamiokande 
detector mostly exist, it does not hold completely. 

Thus, it could be surely concluded from Fig. [11] and Fig. [12] that the SK assumption 
on the direction does not hold as a good estimator for the determination of the directions 
of the incident neutrinos even if statistically. 

Finally, we examine the relation between the zenith angle distribution for upward 
z/ M and and the corresponding zenith angle distribution for // and Jl in the case of no 
oscillation. By perfoming the procedure described in Appendix C, we obtain a pair of 
(cos Op, En) from a sampling of (cos 6^, E v ). In Figure 13, we show the relation between 
the upward incident neutrino zenith angle distribution and the emitted muon ones thus 
obtained for the null oscillation. If we guarantee the SK assumption on the direction, 
the emitted muon zenith angle distribution is expressed approximately as the incident 
neutrino zenith angle distribution. However, the really simulated muon zenith angle 
distribution is quite different from the incident neutrino zenith angle distribution. This 
is the reason why the SK assumption on the direction does not hold even if statistically. 
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Figure 12. Correlation diagrams between cos0„ and cos# p for different muon energy 
ranges. 




Figure 13. The relation between the zenith angle distribution of the incident neutrino 
and corresponding ones of the emitted lepton 



Furthermore, it should be noticed from the figure that the existence of the downward 
muons from the upward neutrinos could not be neglected which is enhanced by the 
azimuthal angles in addition to the backscattering. If the neutrino oscillation does not 
exist, such downward muons do not bring about problems, because in this case the zenith 
angle distribution of the downward neutrino is symmetrical to that of upward neutrino. 
Then, the "inverted leptons" cancel out totally. However, in the case of the presence of 
the neutrino oscillation, such cancellation does not occur. Namely, the determination of 
the direction of the incident neutrinos estimated by the emitted muon should be more 
carefully treated. 
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6. Conclusions 

We have shown the invalidity of the SK assumption on the direction for the analysis of 
the Single Ring \i Events among Fully Contained Events which are solid and free from 
the interpretation on their interaction types. Super-Kamiokande Collaboration have 
examined all possible neutrino events for the neutrino oscillation problem, say, Sub- 
GeV e-like Events, Multi-GeV e-like Events, Sub-GeV /i-like Events, Multi-GeV /x-like 
Events, multi-ring Sub GeV /x-like Events, multi-ring Multi-GeV Events, PC Events, 
Upward Stopping Muon Events and Upward Through Going Muon Events (see, p2]), and 
all possible neutrino events provide the same neutrino oscillation parameters with the 
same precision accidentally. However, different types of neutrino events have different 
structures which are accompanied by different experimental uncertainties. Therefore, 
we could not readily believe such unified conclusions, taking into consideration the 
different experimental qualities in the different types of events. The most clear cut 
events among the SK events are one ring electron like events and muon-like events 
among Fully Contained Events. In these events, apart from the numerical uncertainties, 
one could essentially discriminate electron from muon (see footnote 1), and we decide 
the directions of the leptons as well as their energies, because we know the interaction 
points as well as their end points in the detector. Therefore, it is more desirable that 
one analyze the single-ring electron-like events and single-ring muon-like events among 
Fully Contained Events exclusively not being tempted by the increase of statistics. If one 
really finds the solid evidence for the neutrino oscillation in the most clear cut events, 
such as Single Ring fi Events, then one could find some corroboration for the neutrino 
oscillation even in the worst quality of experimental data, such as Upward Stopping 
Muon Events and Upward Through Going Muon Events. 

In the subsequent paper (Part 2), we apply our present result to the analysis of the 
L/E distribution for Single Ring \i Events due to QEL among Fully Contained Events 
which may lead the direct observation of the neutrino oscillation, if really exists, due to 
the simplicity of the events concerned. Also, we are now under preparation the third 
paper (Part 3) in which we examine the zenith angle distribution for the same type of 
the events mentioned above. 
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In the following Appendices we give the concrete Monte Carlo Simulations, namely, 
the details of our Time Sequential Simulation. 



Appendix A. Appendix: Monte Carlo Procedure for the Decision of 
Emitted Energies of the Leptons and Their Direction Cosines 

Here, we give the Monte Carlo Simulation procedure for obtaining the energy and its 
direction cosines, (l r ,m r ,n r ), of the emitted lepton in QEL for a given energy and its 
direction cosines, (l,m,n), of the incident neutrino. 

The relation among Q 2 , E u+j? , the energy of the incident neutrino, Eg, the energy 
of the emitted lepton (muon or electron or their anti-particles) and 6 S , the scattering 
angle of the emitted lepton, is given as 

Q 2 = 2E v{p) E m {\ - cosfl s ). (A-l) 

Also, the energy of the emitted lepton is given by 

E m = E v{v) - (A-2) 

Procedure 1 

We decide Q 2 from the probability function for the differential cross section with a given 
E v {p) (Eq. (j2J) in the text) by using the uniform random number, £, between (0,1) in 
the following 



where 



i= f Q P e{I) (E u{p) ,Q 2 )dQ 2 , (A-3) 



Pi{£){Eu(u),Q 2 ) = 

&vi{i)(E v (p),Q 2 ) A-Qmax da e ^(E u ( P ), Q 2y 



d Q 2 dQ ■ 



dQ 2 /J QLn dQ 2 

(A-4) 



From Eq. (A-l), we obtain Q 2 in histograms together with the corresponding theoret- 
ical curve in Fig. IA11 The agreement between the sampling data and the theoretical 
curve is excellent, which shows the validity of the utlized procedure in Eq. (A -3) is right. 

Procedure 2 

We obtain E e ^ from Eq. (A-2) for the given E u ^ and Q 2 thus decided in the Procedure 
1. 



Procedure 3 

We obtain cos# s , cosine of the the scattering angle of the emitted lepton, for E e ^ thus 




Figure Al. The reappearance of the probability function for QEL cross section. 
Histograms are sampling results, while the curves concerned are theoretical ones for 
given incident energies. 



decided in the Procedure 2 from Eq. (A-l) . 



Procedure 4 

We decide 0, the azimuthal angle of the scattering lepton, which is obtained from 



(p = 2tt£ . (A-5) 
Here, £ is a uniform random number of the range (0, 1). 

As explained schematically in the text (see Fig. [3] in the text), we must take account of 
the effect due to the azimuthal angle in the QEL to obtain the zenith angle distribu- 
tion of both Fully Contained Events and Partially Contained Events correctly. 

Procedure 5 

The relation between direction cosines of the incident neutrinos, (£ u (u), m u(u), n u(u)), and 
those of the corresponding emitted lepton, (£ T , m r , n r ), for a certain 9 S and is given as 



hi 



rn 



( ir \ 








V n * ) 





-v{v) 



mn 



o 



m 



n u (p) J 



' sin# s cos0 \ 

sin# s sin0 
V cos6» s , J 



(A-6) 



where rwp) = cos9 u ^, and n r = cos0£. Here, 9e is the zenith angle of the emitted lepton. 



The Monte Carlo procedure for the determination of 9g of the emitted lepton for 
the parent (anti-)neutrino with given 9 v (p) and E v {p) involves the following steps: 
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Figure A2. The relation between the direction cosine of the incident neutrino and 
that of the emitted charged lepton. 



We obtain (£ r , m r , n r ) by using Eq. (1A-6I) . The n r is the cosine of the zenith angle 
of the emitted lepton which should be contrasted to n u , that of the incident neutrino. 
Repeating the procedures 1 to 5 just mentioned above, we obtain the zenith angle dis- 
tribution of the emitted leptons for a given zenth angle of the incident neutrino with a 
definite energy. 

In the SK analysis, instead of Eq. flA-6j) . they assume 
n r = n^( P ) uniquely for E^cp) > 400 MeV. 



Appendix B. Appendix: Monte Carlo Procedure to Obtain the Zenith 
Angle of the Emitted Lepton for a Given Zentith Angle of the Incident 
Neutrino 

The present simulation procedure for a given zenith angle of the incident neutrino starts 
from the atmospheric neutrino spectrum at the opposite site of the Earth to the SK 
detector. We define, Ni nt (E u ^, t, cos9 v ^\), the interaction neutrino spectrum at the 
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depth t from the SK detector in the following way 

N int (E u{p) ,t,cos9 u{D) ) = N sp (E„ (p) ,cos9 u{p) )x 

( dt \ / dt \ 

x • • • x II — 



^l(E v (v), ti, pi) J V ^n\E v (v)i tni P; 



(b-i; 



Here, N S p(E u ^, cos^p)) is the atmospheric (anti-) neutrino spectrum for the zenith 
angle at the opposite surface of the Earth. 

Here Xi(E u ^, U, pi) denotes the mean free path of the neutrino (anti neutrino) with 
the energy E v ^ due to QEL at the distance, t iy from the opposite surface of the Earth 
inside whose density is p^. 

The procedures of the Monte Carlo Simulation for the incident neutrino (anti neu- 
trino) with a given energy, £v(p), whose incident direction is expressde by (/,m, n) are 
as follows. 



Procedure A 

For the given zenith angle of the incident neutrino, V ^, we formulate, 
Np TO (E u ( j? ),t,cos9 u ( 9 ))dE u (p} , the production function for the neutrino flux to produce 
leptons at the Kamioka site in the following 

N pro (E v ( P -),t, cos 9 v { y v))dE vly9 ) 

— o~i(i)(Eis(D))N in t(E l/ ( D ),t,cos9 u ( D ))dE u (p), (B-2) 

where 



C max 



^)(^(*)) = (> dQ V d Q- (B-3) 

Each differential cross section above is given in Eq. (T5]) in the text. 

Utilizing, £, the uniform random number between (0,1), we determine E v (p\, the energy 

of the incident neutrino in the following sampling procedure 



6 = / " M P d {E u{D)) t )C os9 u{D) )dE u{D)l (B-4) 

^v(is) ,min 



where 



Pd(E„(o), t, cos9 l/ ( 9 ))dE l/ ( 9 ) 

N P ro(E v ( D ),t, cos9 u ( 9 ))dE v ( V ) 



f(f),max 

N wo {E v (p) , t , cos 9 v{ p) ) dE v (p) 



,min 



In our Monte Carlo procedure, 
the reproduction of, P^{E v i^ t, cos6 v (p))dE v fp\, the normalized differential neutrino in- 
teraction probability function, is confirmed in the same way as in Eq. (A-4). 
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Procedure B 

For the (anti-)neutrino concerned with the energy of E V ( D \, we sample Q 2 utlizing £ 3 , 
the uniform random number between (0,1). The Procedure B is exactly the same as in 
the Procedure 1 in the Appendix A. 

Procedure C 

We decide, 9 S , the scattering angle of the emitted lepton for given E v r p \ and Q 2 . The 
procedure C is exactly the same as in the combination of Procedures 2 and 3 in the 
Appendix A. 

Procedure D 

We randomly sample the azimuthal angle of the charged lepton concerned. The Proce- 
dure D is exactly the same as in the Procedure 4 in the Appendix A. 

Procedure E 

We decide the direction cosine of the charged lepton concerned. The Procedure E is 
exactly the same as in the Procedure 5 in the Appendix A. 

We repeat Procedures A to E until we reach the desired trial number. 



Appendix C. Appendix: Correlation between the Zenith Angles of the 
Incident Neutrinos and Those of the Emitted Leptons 

Procedure A 

By using, N pro (E u ^,t, cos 9 u ^))dE u ^, which is defined in Eq. (IB-2j) . 
we define the spectrum for cos 6 v (p) in the following. 



J(cos6^(p))d(cos6^( P )) = 

d(cos0„(ij)) / N pro (E u{p) ,t, cos9 u{p) )dE u{p) . (C-l) 

By using Eq. ( 1C-2I) and £, a sampled uniform random number between (0,1), then we 
could determine cos 9 v (p\ from the following equation 

/■COS B v {n\ 

t = / P n (cos9 l/(p) )d(cos9 l/(p) ), (C-2) 



where 



P„(cos6 l I/ (p)) = I(cos9 v (p)) I I(cos9 v (p))<\{cos9 v{ p)). (C-3) 
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Procedure B 

For the sampled d(cos9 u ^) in the Procedure A, we sample E v t D ) from Eq. flC-41) by using 
£, the uniform randum number between (0,1) 



P wo {E v{p ), cos w (p))dE„( P ), (C4) 



where 



P pr . (E I/ ( P ),t,cos6 l I/ (p))dE ! ,(p) = 

N pro (Ey(„) , t , COS fl^p) )dE„(jj) 



'^(fjjinax 



^(y},min 



NproiE^p)^, cos I/ ( j7 ))d£' i/ (p) 



(C-5) 



Procedure C 

For the sampled -Ev(p) in the Procedure B, we sample E^qa from Eqs. (IA-2[) and (IA-3|) . 
For the sampled E u ^ and E^q^, we determine cos# s , the scattering angle of the muon 
uniquely from Eq. (1A-1I) . 

Procedure D 

We determine, 0, the azimuthal angle of the scattering lepton from Eq. (IA-5[) by using 
£, an uniform randum number between (0,1). 

Procedure E 

We obtain cos from Eq. (IA-61) . As the result, we obtain a pair of (cos 9 U ^, cos 6^)) 
through Procedures A to E. Repeating the Procedures A to E, we finally the correlation 
between the zenith angle of the incident neutrino and that of the emitted muon. 
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